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Frequency dependence of the large frequency separation of 
solar-like oscillators: Influence of the Helium 
second-ionization zone 
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ABSTRACT 

The large frequency separation (Ai/) between modes of the same degree and consec- 
utive orders in a star is approximately proportional to the square root of its mean 
density. To determine Av as accurately as possible a mean large frequency separation 
((Az/)) computed over several orders is often used. It is, however, known that Av 
varies with frequency in a second order effect. From observations it has been shown 
that this frequency dependence is more important for main-sequence stars than it is 
for red-giant stars. Here we use YREC models to verify and explain this observational 
result. We find that for stars with R > 8 R© the effect of the Helium second ionisation 
zone is relatively small. For these stars the deep location of the He II zone induces a 
frequency modulation covering only a few Av, while the amplitude of the modulation 
is low due to the relatively weak and extended He II layer, causing a shallow wide 
depression in Ti. For less evolved stars the He II zone is located closer to the surface, 
and it is more confined, i.e. a deep narrow depression in I\. This causes frequency 
modulations with relatively high amplitudes covering up to about 20Av, inducing a 
relatively large frequency modulation. Additionally, we find that for less evolved stars 
the He II zone is stronger and more localised for more massive stars and for stars with 
low metallicities further increasing the amplitude of the frequency modulation. 
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1 INTRODUCTION 

Stellar oscillations can be used to determine the internal 
structures of stars. Solar-type stars, subgiants and red-giant 
stars exhibit solar-like oscillations, i.e., oscillations stochas- 
tically excited in the turbulent outer layer. These oscillations 
generally appear in regular patterns in frequency following 
to reasonab l e ap proximation the asymptotic relation derived 
bv lTassoull (ll980D : 



u n>l « Av{n + 1/2 + e) - £(£ + 1)D , 



(1) 



with v„ t i the frequency of an oscillation mode with radial 
order n and degree I and Av the large frequency separation 
between modes of the same degree and consecutive orders. 
Do is most sensitive to deeper layers in the star and e to the 
surface layers. 

The large frequency separation is inversely propor- 
tional to the sound travel time through the star and it 
can be shown that this is directly proport ional to the 
square root of the mean density of the star (jUlrichl Fl986l : 
iKieldsen fe Bedding! [l995h . To a first order approximation 



the large frequency separation is constant over the observed 
frequency range and often a mean large frequency separa- 
tion ((Av)) is computed. However, for the Sun and other 
main-sequence stars it has been shown that the frequency 
range over which (Av) is computed has a significant in- 
fluence on the resulting value. This frequency dependence 
can be attributed to acoustic glitches, i.e., sudden inter- 
nal property changes at the base of the convection zone 
and at the Helium second-ionization zone. Also other more 
slowly varying underlying variations, such as changing con- 
ditions close to the core, can cause departures from uni- 
form Av. Theory and methods considering acoustic glitches 
were first explored in the conte xt of studies o f solar oscil- 
lations, as initially proposed in iGoughl (Il990|). and devel- 
oped fu rther by e.g. Basu fc Antial lll994l), llVIonteiro et al.l 
1 19941). Perez Hernandez fc Christensen-Dalsgaardl (|l99^ )T 
iBasu fc Mandell (|2004n .' Different methods of analysis for 
extracting data on the glitches h av e been developed 
(e.g. see iMazumdar fc Antial l2001a|lbl: iBallot et al l |2004| ; 



IVerner et al.ll2006l : lHoudek fc Goughll2007l : lRoxburghll2010l ). 



© 2011 RAS 



2 S. Hekker et al. 




ants using mod els constructed using Y REC, the Yale stellar 
evolution code (|Demarque et al]|200St ). 



Figure 1. YREC evolutionary tracks in a Hertzsprung- Russell di- 
agram. Tracks with approximately solar metallicity are indicated 
with diamonds, while the solid lines show the tracks with [Fe/H] 
= +0.3 and the dashed lines indicate the tracks with [Fe/H] = 
-0.5 for models with M = 1.0 M (black) and 2.4 M (red). 
Models with M = 1.0 M , [Fe/H] = 0.0 including either dif- 
fusion, Krishna Swamy T-r relation or a = 1.75 H p (see text 
for more details) are indicated with the black, mint and purple 
dashed-dotted lines, respectively. These lines nearly coincide with 
the solar metallicity M = 1.0 M models indicated with the black 
diamonds. The black cross indicates the model used in Fig. [2] 
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Figure 2. Au computed using individual frequencies of I = 
modes as a function of radial order (bottom axis) or frequency 
(top axis) for the 1.2 M star indicated with the black cross in 
Fig.[l] The intervals used to compute (Au) are indicated with the 
red and green vertical lines. 



in some cases for application to other stars. Recently, the 
Hell ionization zon e glitch signal has been measured in th e 
red giant H R7349 dCarrier et al.ll20ld: iMiglio et al.ll20ici l. 
Additionally, iMazumdar fc Michel! I|2010l ) have claimed de- 
tection of acoustic glitches, from the Hell ionization zone 
and the base of the convective envelope, in solar-like oscil- 
lations of the F-type star HD49933. 

In contrast to the frequency dependence of {Au) for the 
Sun and main-sequence star s, first results from the Kepler 
mission l|Borucki et alj|20ich showed that for red giants the 
value of (Au) is much less se nsitive to the freque ncy range 
over which it was computed (| Hekker et al.|[201ll ). Here we 
investigate the reason for the difference in sensitivity of (Au) 
to the frequency range for main-sequence stars and red gi- 



be found 
OPAL opacities 



2 YREC MODELS 

A full overview o f the YREC code can 
in Dem arquc c t al. (120081 ). We use 
ijlglesias fc Rogers! 1996) supple mented with low tem pera- 
ture (logT < 4.1) opacities of iFerguson et al.l |2005l) and 
the OPAL equation of state ^Rogers fc Nayfonovll2002l). All 
nuclea r reaction rates are obtained from lAdelberger et all 
ljl998t) . except for that of the 14 N(p,-y ) 15 Q r eaction, for 
which we use the rate of iFormicola et all (|2004T ). 

In this work we use a sequence of models of solar metal- 
licity with masses between 1.0 and 2.4 M . The sequences 
range from just after ZAMS to near the tip of the giant 
branch for stars of masses up to 2.0 M and up to the he- 
lium burning stage for the higher mass models. Additionally 
we also use 1.0 and 2.4 M models with [Fe/H]=+0.3 and 
—0.5 to investigate the effects of metallicity. The masses and 
metallicities are chosen to bracket the masses and metal- 
licities of stars expected to be observed with Kepler. The 
models were constructed without diffusion and gravitational 
settling of helium and heavy elements. However, we do ex- 
am ine the role of diffu sion and gravitational settling using 
the lThoul et al.l l|l994h prescription for a sequence of 1.0 M 
models. All models described here were constructed to have 
the Eddington T-r relation in their atmosphere. To study 
the effect of atmospheric structure we use a sequence of TO 
M m odels constructed with the Krishna Swamy T-r re- 
lation (|Krishna Swamv| [l966). Additionally, to examine the 
effect of the mixing-length parameter, a, we have made a 
sequence of 1.0 M models with a = 1.75 H p instead of 
the solar value of a = 1.826 H p used in the other models. 
Evolutionary tracks of the models are shown in Fig. [T] 

Note that we consider solar-like oscillations in models 
of main-sequence stars up to masses of 2.4 M , although 
solar-like oscillations are generally not observed in main- 
sequence stars with masses above roughly 1.5 M . Further- 
more, we investigate models with radii less than 20 R and 
wi th Au > 2 ^H z com patible with the observational results 



wi tn l\v > I jitV L compa 
bv lHekker etaT] (|201ll ). 



DEPENDENCE OF {Au) 
RANGE 



ON FREQUENCY 



To investigate the frequency dependence of {Au) from the 
models, we derived {Au) from linear fits of the frequen- 
cies of radial modes in the intervals f max ± 2(A!/) sca ii ng 
and ^ m ax ± 4(Ai/) sca ii ng as a function of their radial orders 
(see Fig. [2] and Eq. [I}. ^ max is the frequency of maximum 
oscillation power estimated from t he mass, radius an d ef- 
fective temperature of the models l|Brown et al.|[l99ll ) and 
(A^) sca ii ng is the mean large separation obtained fr om the 
mass and radius of the star, usi ng scaling relations (|Ulrichl 
1 19861 ; iKieldsen fc Bedding! Il99sh . We also investigated the 
median Au values in the respective frequency intervals to 
investigate the symmetry of Au in frequency. This showed 
that the behaviour of Au with frequency is antisymmetric 
about Vmax and that a linear fit can be used in a first order 
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Figure 3. Left: Ratio of (Av) values computed from a linear fit to the frequencies versus the radial order of the modes in the intervals 
fmax i 2Af and v max ± 4Af as a function of stellar radius for different masses (top) and different metallicities or diffusion (bottom). 
Right: Same as right panels but now as a as a function of (Av). The dashed black line indicates unity. 



no approximation. We note that in those cases where it is not 
possible to extract individual radial-mode frequencies from 
real data, techniques are usually applied to the full oscil- 
lation power spectrum to estimate (Av), in which case the 
estimated value is also affected by non-radial modes. 

115 The ratio of the (Av) values computed over each fre- 

quency range gives an indication of the sensitivity of (Av) to 
the frequency range. Fig. [3] shows this ratio as a function of 
stellar radius. The scatter in the ratio of (Av) computed over 
different intervals symmetric around v m „ decreases with in- 

120 creasing stellar radii and that the ratios converge to 1.0. 

From the ratios of (Az/) we find that the scatter at low 
radii (< 8 R0) reduces for masses below roughly 1.5 Mq. 
On the other hand no significant change in the scatter is 
present due to different metallicities or due to the inclusion 

125 of diffusion, different atmospheres or mixing-length param- 
eters in the models (see bottom panels of Fig. [2]). Note that 
the relative change in the values of (Ai/) for stars with radii 
< 8 R,0 as a function of the frequency range is of the same 
order as the precision as with which (Az/) can be determined 

130 from current state-of-the-art data. 

iHekker et ail (|201lh suggested two possible reasons for 
the decreased influence of the frequency range on (Ai/) for 
red giants: 1) the trend in Av over a typical frequency range 
is approximately linear, and/or 2) Af changes relatively 

135 slowly with frequency. Other reasons might be that in red gi- 
ants fewer modes are observed, and these modes have lower 
radial orders than in less evolved stars. We test these sug- 
gestions by fitting a linear polynomial through the values of 



Av over a range ;-'max±4(A^) sca ii ng and investigate the slope 
of the fit and the standard deviation of the values around 140 
the fit. These results are shown in Fig. 0] and show that 
the slope or linear trend of the variation in Av as a func- 
tion of frequency for stars with lower Av, i.e., stars with 
larger radii, is larger than for stars with higher values of 
Av, with little sensitivity to the mass, metallicity, diffusion, 145 
atmosphere and mixing-length parameter. For higher mass 
models there is some variation visible in the trend, which 
coincides with the 'hook' in the H-R diagram (Fig. [!}, i.e., 
the contraction during the turn-off of the main sequence. 
The scatter around the fit is however much lower for stars 150 
with Av below ~10 /iHz than for stars with larger values 
of Av. Thus we indeed conclude that stars with lower Av 
(< 10 /^Hz), i.e., larger radii (> 8 R©), have a predominantly 
linear dependence on frequency, while for stars with larger 
Av (> 10 /iHz), i.e., smaller radii (< 8 R©), the trend is 155 
shallower and the scatter in Av increases. The amount of 
scatter does not seem to depend critically on the inclusion 
of diffusion, different atmosphere or mixing-length param- 
eter in the models, while it increases with mass and with 
decreasing metallicity. 160 



4 DISCUSSION 

Acoustic glitches, i.e., regions of sharp-structure variation 
in the stellar interior, are known to cause a modulation of 
Av with frequency. For the Sun and other main-sequence 
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Figure 4. Left: Linear trend in Av over the range i/ ma x±4Ac for models with different masses (top) and different metallicity, diffusion, 
atmosphere or mixing-length parameter (bottom) as a function of Ac derived from scaling relations. Right: Standard deviation of Av 
after subtracting the linear fit over the range c max ± 4Av. 



stars both the Helium second-ionization zone (He II zone) 
and the base of the convection zone have been shown 
to contribute to the modulation of the freque ncies (e.g. 
iMonteiro fc Thompson! 120051 ; IVerner et al.ll200d and refer- 
ences in the introduction). Generally, the signal of the base 
of the convection zone is weaker than the signal from the He 
II zone, due to its location deeper in the star. Additionally, 
iMiglio et alJ {2010) already showed that for a red giant only 
the He II zone caused a measurable modulation. 

The He II zone causes a local depression in the first 
adiabatic exponent 71 = (d lnp/dln p) s , with p pressure, 
p density and s specific entropy. This glitch modulates the 
frequencies in a sinusoidal manner with a so-called 'period' 
inversely proportional to the acoustic depth (r), with 




in which r indicates the radius of the glitch, R the ra- 
dius of the star and c sound speed. The 'period' due to 
the bottom of the convection zone, located at the posi- 
tion where VT becomes larger than (VT) a a, has been com- 
puted similarly with respect to the closest acoustic boundary 
i|Mazumdar fc Antidl200l"bl) . 

The period due to the He II zone depends on the loca- 
tion of the depression in 71, i.e., a depression close to the 
surface will cause a lower value for r than a depression lo- 
cated deeper in the star. We computed the expected period 
(t -1 ) of the frequency modulation from the location of the 
He II zone for all our models and divided it by (A^) sca ii ng to 



investigate how many orders one modulation period covers. 
This is shown in Fig. [5] For models with Av < 10 ^iHz the 
modulation period is about 5Av and not significantly depen- 
dent on mass, metallicity, diffusion, atmosphere or mixing- 
length parameter, while the period and the spread in periods 
increase for models with Av > 10 /iHz. 

To determine the dependence of (Ai/) on the frequency 
range we added/removed 2Av on either side of the frequency 
interval (Section 3). This is about one full period for models 
with Av < 10 fiRz and no significant net effect in the mean 
value of Av has been observed. For stars with Av > 10 ^tHz 
the period depends on mass and is in some cases signifi- 
cantly longer than hAv implying that only part of a period 
is added/removed for the different frequency ranges. This 
does have an effect on the resulting mean values of Av. From 
this we can conclude that for larger stars (R > 8 R©) the 
smooth stellar structure changes reflected in the linear trend 
of Av are dominant over the effect from the He II acoustic 
glitch, while for smaller stars (R < 8 R©) the effect from 
the He II acoustic glitch is dominant. 

In addition to the period, we also investigated the am- 
plitude of the variations in Av, parametrized by the stan- 
dard deviation (see right hand side of Fig. [4}. The amplitude 
of the modulation is larger for stars with Av > 10 /iHz than 
for stars with Av < 10 ^iHz and depends on mass and metal- 
licity (see Fig. [4}. We expect the amplitude to depend on the 
shape of the glitch. We parametrize the shape of the glitch 
190 by the ratio of the depth of the He II depression in 71, i.e. 
220 its strength, to the width (in units of acoustic radius) of the 
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Figure 5. Period of the frequency modulation expressed in units 
of Av as a function of Av. 
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Figure 7. For M = 1.2 Mq models the standard deviation of the 
variation in Av (asterisks) and the depth of the He II depression 
in 7i divided by the width of the depression expressed in acoustic 
radius (diamonds) as a function of Av. 
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Figure 6. The depth of the He II depression in 71 divided by the 
width of the depression expressed in acoustic radius as a function 
of Av. 



depression, i.e., the narrower the width for a given depth 
the bigger the effect on the frequencies. The ratios for all 
models are shown in Fig. [6] as a function of Av. 

The changes in the shape of the He II depression in 71 

225 at a given Av are clearly reflected in the changes in the stan- 
dard deviation of the variation in Av at that Av (see Fig. [7]) . 
The general trends in the curves for Av and the depth to 
width ratio of the He II zone are the same, although there 
are differences in the amplitude of the variations. Therefore 

230 we conclude that the shape of the He II depression in 71, i.e. 
the ratio of the strength over the width of the He II zone, 
is a significant cause of the difference in the amplitudes of 
the variation in Av which vary with mass and metallicity. 
However, there must be additional effects that play a role. 

235 This could be the base of the convection zone, which induces 
a small modulation with a period that is of the same order 
as the period of the He II zone. 

We note that the increasing value of Av at which the 
amplitude of the variation in Av increases is a function of 
mass and metallicity (Fig. Q and that this effect is consis- 240 
tently present in both the period and shape of the He II zone 
(Figs.[5]&[6]). 



5 CONCLUSIONS 

In th is work we confirm the statement by iHekker et al.l 
l|201ll ) that the value of (Av) computed for red giants is 245 
not as sensitive to the frequency range over which it is com- 
puted as is the case for the Sun or other main-sequence 
stars. The analysis of frequencies from YREC models of 
stars with different masses and metallicities and with the 
inclusion of diffusion, a different atmosphere and different 250 
mixing-length parameter clearly show that for stars with 
lower Av (< 10 piHz), i.e., larger radii, the dependence of 
Av follows a linear trend with relatively low scatter. For 
stars with higher Av (> 10 A*Hz), i.e., smaller radii, the 
linear trend diminishes, while the scatter increases, with a 255 
peak in the scatter at about 100 pHz. The magnitude of this 
peak in the scatter increases with mass and with decreasing 
metallicity, but does not seem to be sensitive to the inclu- 
sion of diffusion, a different atmosphere or mixing-length 
parameter in the models. 260 

The lower sensitivity of (Av) on frequency is due to 
deeper location of a weaker and wider He II zone for ex- 
tended evolved stars. For less evolved stars the shallower 
location of a stronger and thinner He II zone causes an in- 
creased sensitivity of (Av) on frequency. The strength and 265 
width of the He II zone in less evolved stars changes as a 
function of mass and metallicity. 



ACKNOWLEDGEMENTS 

SH acknowledges financial support from the Netherlands Or- 
ganisation for Scientific Research (NWO). YE and WJC ac- 270 
knowledge financial support of UK STFC. SH thanks M. 
Mooij for useful discussions. 



REFERENCES 

Adelberger E. G., et al. 1998, Reviews of Modern Physics, 
70, 1265 

Ballot J., Turck-Chieze S., Garcia R. A., 2004, A&A, 423, 
1051 

Basu S., Antia H. M., 1994, MNRAS, 269, 1137 



© 2011 RAS, MNRAS 000, \JM 



6 S. Hekker et al. 



Basu S., Mandel A., 2004, ApJ, 617, L155 

280 Borucki W. J., et al. 2010, Science, 327, 977 
Brown T. M., et al. 1991, ApJ, 368, 599 
Carrier F., et al. 2010, A&A, 509, A73 
Demarque P., et al. 2008, Ap&SS, 316, 31 
Ferguson J. W., et al. 2005, ApJ, 623, 585 

285 Formicola A., et al. 2004, Physics Letters B, 591, 61 
Gough D. O., 1990, Nature, 347, 346 
Hekker S., et al. 2011, A&A, 525, A131 
Houdek G., Gough D. O., 2007, MNRAS, 375, 861 
Iglesias C. A., Rogers F. J., 1996, ApJ, 464, 943 

290 Kjeldsen H., Bedding T. R., 1995, A&A, 293, 87 
Krishna Swamy K. S., 1966, ApJ, 145, 174 
Mazumdar A., Antia H. M., 2001a, A&A, 368, L8 
Mazumdar A., Antia H. M., 2001b, A&A, 377, 192 
Mazumdar A., Michel E., 2010, ArXiv: 1004.2739 
Miglio A., et al. 2010, A&A, 520, L6 295 
Monteiro M. J. P. F. G., Christensen-Dalsgaard J., Thomp- 
son M. J., 1994, A&A, 283, 247 
Monteiro M. J. P. F. G., Thompson M. J., 2005, MNRAS, 
361, 1187 

Perez Hernandez F., Christensen-Dalsgaard J., 1998, MN- 300 

RAS, 295, 344 
Rogers F. J., Nayfonov A., 2002, ApJ, 576, 1064 
Roxburgh I. W., 2010, Ap&SS, 328, 3 
Tassoul M., 1980, ApJS, 43, 469 

Thoul A. A., Bahcall J. N., Loeb A., 1994, ApJ, 421, 828 305 
Ulrich R. K., 1986, ApJ, 306, L37 

Verner G. A., Chaplin W. J., Elsworth Y., 2006, ApJ, 638, 
440 



© 2011 RAS, MNRAS 000, [TUB] 



